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We find that the chromomagnetic instability existing in neutral two- flavor color superconductivity 
at moderate densities is removed by the formation of an inhomogeneous condensate of charged gluons 
and the corresponding induction of a magnetic field. It is shown that this inhomogeneous ground 
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in a color superconductor at moderate densities can be of interest for the astrophysics of compact 
stellar objects exhibiting strong magnetic fields as magnetars. 
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I. INTRODUCTION 

It is quite plausible that color superconductivity (CS) 
is realized in the cold and dense core of compact stellar 
objects. For matter densities high enough to guarantee 
the quark deconfinement, and sufficiently low to consider 
the decoupling of the strange quark mass M s , the quarks 
can condense to form a two-flavor color superconductor 
(2SC) dl|. At higher densities, for which M s cannot be re- 
garded as extremely large, the three-flavor phases gCFL 
[H and CFL Q will appear respectively with increasing 
densities. Matter inside compact stars should be neu- 
tral and remain in (3 equilibrium. When these conditions 
along with M s are taken into account in the dynamics 
of the gluons in the gCFL phase (or for that matter in 
the 2SC phase if the s-quark is decoupled), some gluon 
modes become tachyonic 0, 0] , indicating that the sys- 
tem ground state should be restructured. 

The quest to find the stable ground state at the realis- 
tic, moderate densities of compact stars has led to several 
interesting proposals. Some of the most promising pos- 
sibilities are a modified CFL-phase with a condensate 
of kaons that requires a space-dependent condensate 
to remove the instability; a LOFF phase on which the 
quarks pair with nonzero total momentum Q ; and a ho- 
mogenous gluon condensate phase At present it is 
not clear if any of these proposals is the final solution to 
the problem. 

In this paper we will show that there is one more pos- 
sible solution: an inhomogeneous condensate of charged 
gluons that induces a rotated magnetic field. Our so- 
lution is found in the Meissner unstable region of the 
so-called gapped 2SC phase [§], but we expect that it 
could also be realized in the three-flavor theory. It has 
the peculiarity of preserving the electromagnetic gauge 
invariance U em (l) of the color superconductor. As will 
be shown in our derivations, this inhomogeneous phase 
is energetically favored over the one with homogeneous 
gluon condensate and no induced magnetic field. We will 
also discuss its possible relation to a LOFF-type phase. 

In a previous paper (To| . we investigated the forma- 
tion of an inhomogeneous gluon condensate in a color 
superconductor. However, apart from the fact that in the 



current work we a re g oing to consider a two-flavor the- 
ory, while in Ref. [10( we studied a three-flavor system, 
there is a fundamental difference between the genesis of 
the condensation phenomenon that we will discuss in the 
present paper and the one found in [10]. In Ref. 
the condensation of gluons was the response of the color 
superconductor to an externally applied and sufficiently 
strong magnetic field. In that case, if no magnetic field 
is present, no gluon condensate is formed. Therefore, the 
gluon condensation in [l(| was a magnetic-field-induced 
phenomenon. On the other hand, in this work, no ex- 
ternal magnetic field is present to begin with. The con- 
densation phenomenon here is connected to a Meissner 
instability triggered at moderate densities by a pairing 
stress associated with the neutrality and /3-equilibrium 
constraints. Such instability occurs in the gapped 2SC 
color superconductor at moderate densities when the ef- 
fective Meissner mass of the charged gluons becomes 
tachyonic. As we shall show below, in this situation an 
inhomogeneous condensate of charged gluons emerges to 
remove the chromomagnetic instability created by the 
pairing mismatch. 

In the case of the gluon condensation produced by 
a sufficiently strong magnetic field [hJ, the anomalous 
magnetic moment interaction of the spin-1 charged fields 
with the external magnetic field led to a phenomenon of 
anti-screening, so that the back reaction of the conden- 
sate on the applied field was such that it boosted the in- 
medium magnetic field strength. This anti-screening is a 
well-known phenomenon that occurs in all the instances 
of charged spin-1 field condensation in the presence of 
a strong magnetic field [ll], [ijj ■ Amazingly, as we will 
show in this paper, the system predisposition to boost 
an in-medium magnetic field through the magnetic mo- 
ment of the spin-1 charged condensate works in favor of 
the spontaneous generation of a rotated magnetic field 
in the unstable region of the gapped 2SC when no ex- 
ternal magnetic field is present. Indeed, in this case the 
currents associated to the inhomogeneous condensate of 
charged gluons induce a non-zero rotated magnetic field 
in the system. Hence, the magnetic field here is a con- 
sequence of the gluon condensation, while in 10] it was 
the other way around. The spontaneous induction of a 
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rotated magnetic field in CS systems that have pairings 
with mismatched Fermi surfaces is a new kind of phe- 
nomenon that can serve to generate magnetic fields in 
stellar compact objects as magnetars. 

The plan of the paper is as follows. In Sec|TT]we present 
the effective theory of the gluonic phase with electromag- 
netic interactions in the neutral two- flavor dense QCD. 
The arising of an inhomogeneous gluon condensate to 
remove the chromomagnetic instability caused by the 
Fermi-surface pairing mismatch in the 2SC color super- 
conductor at moderate densities, as well as the corre- 
sponding induction of a rotated magnetic field are inves- 
tigated in Sec lIIII Finally, in Sec. IIVI we discuss the 
significance of the reported findings for the astrophysics 
of highly magnetized compact objects. 



II. NEUTRAL TWO-FLAVOR DENSE QCD 
WITH ELECTROMAGNETIC INTERACTIONS 

In the gapped 2SG phase the solution of the neutral- 
ity conditions dflo/dfii — 0, with fa — /Lt e ,/X8,/i3, an d 
gap equation dCi a /dA — 0, for the effective potential 
fio in the mean-field approximation, led to /Z3 = 0, and 
nonzero values of fi e , and /j,g, satisfying /ig <C He < M 
for a wide range of parameters @ . Here [i is the quark 
chemical potential, /i e the electric chemical potential, and 
the "chemical potentials" /ig and /13 are strictly speak- 
ing condensates of the time components of gauge fields, 

Ms = (\/3 5 /2)(G^ 8) ) and ^3 = (.9/2)(G^ 3) ). The nonzero 
values of the chemical potentials produce a mismatch 
between the Fermi spheres of the quark Cooper pairs, 

The gapped 2SC turned out to be unstable once 
the gauge fields , G (2) , , K„ K\ , & A^} were 
taken into consideration. As shown in Ref. j4j, the glu- 

ons G^ 1 ' 2 ' 3 - 1 are massless, the in-medium 8 tft, -gluon G^ = 
sin OA^ + cos 8 G^ has positive Meissner square mass, and 
the X-gluon doublet JfT = (G^ - iGf , g£ 6) - iG$P ) 
has Meissner square mass that becomes imaginary for 
A > 5(i > A/y2, signalizing the onset of an unstable 
ground state. The mass of the in-medium (rotated) elec- 
tromagnetic field A^ — cos 8 A^ — sin#G^ is zero, which 
is consistent with the remaining unbroken t/(l) em group. 
Note that even though the gluons are neutral with re- 
spect to the regular electric charge, in the 2SC phase the 
gluon complex doublet has nonzero rotated charge 
Q = Q — ^= T§ of magnitude q — e/2, with e = e cos 8, Tg 
the 8th color generator, and Q the conventional electric 
charge operator [l3j |. 

In what follows, we will find a stable ground state so- 
lution near the critical point Sfj, c = A/v2- Close to that 
point the absolute value of the square magnetic mass 
\m 2 M I is very small, thus allowing analytical calculations. 
The tachyonic modes in the Meissner-unstable gapped 
2SC are associated only to charged gluon fields. It is nat- 



ural to expect then that a new, stable ground state should 
incorporate the condensation of these charged gluons. In 
general we should allow for an inhomogeneous conden- 
sate. Given that this kind of solution may generate ro- 
tated electromagnetic currents, the rotated electromag- 
netism should be included in the general framework of the 
condensation phenomenon. With this aim, we start from 
the mean-field effective action of the 2SG gauged Nambu- 
Jona-Lasinio model depending on the gluon fields, ro- 
tated electromagnetic field, diquark condensates, and 
chemical potentials, 

r e// = -/ d' x {\{f^f + I ( F;j 2 }+r 9 -^^+r g) 

(i) 

where a = 1,...,8, G is the four-fermion coupling, 
(3 the inverse temperature, V the 3-dimensional vol- 
ume, f^ v = d^Au — d v A^ is the conventional elec- 
tromagnetic strength tensor, T g is the gauge-fixing ac- 
tion, and T q = iTrln(5'- 1 + G^ + A^) is the 
fermion contribution obtained after integrating out the 
quark fields in the grand partition function of the 
neutral 2SG phase [Tj]. In T q , S is the quasipar- 
ticle propagator in the Nambu-Gor'kov space; T£ — 
diag^gj^Ta, —gj^TJ) is the Nambu-Gor'kov gluon ver- 
tex, and r M = diag(e^ fl Q, —ej^Q) is the Nambu-Gor'kov 
electromagnetic vertex. In principle, S depends on the di- 
quark condensate, chemical potentials and possible con- 
densates of gauge fields. 

We assume a weak-coupling approximation ( a s < 1). 
Then, expanding the logarithm up to second order in the 
gauge fields, we obtain r, = T - \ J d 4 x[G°n% G b v + 
AfJI^" A v ]. All the quantities contributing to T q : the 
effective action (Tq = ^Tr In S -1 ), as well as the polar- 
ization operators for gluons = ^Tr[T^ST^S]) and 
photons (II^ = i TV [r^ 1 ST" SI) have been found in the 
hard-dense-loop approximation of the neutral 2SC phase 
[J]. The absence of linear terms in the expansion of T q 
is due to the zero trace in the color index, which is a 
consequence of color neutrality, (J c ) = (i.e. zero tad- 
pole), required to preserve the theory non-Abelian gauge 
symmetry [l5j . 

In this derivation we used t'Hooft gauge to avoid mix- 
ing between the gluons and the fluctuations of the di- 
quark condensate. Details of the implementation of 
t'Hooft gauge in the 2SC theory can be found in [Til ]. 
The scalar modes associated to the fluctuations of the 
diquark condensate come from the spontaneous breaking 
of the color gauge group. Therefore, they are unphys- 
ical and only serve to provide the longitudinal degrees 
of freedom of the gluons that acquire a mass due to the 
symmetry breaking. The decoupling of the scalar modes 
from the gluons occurs in the t'Hooft gauge for arbitrary 
values of the gauge parameter A. Later on we will choose 
A = 1 because it simplifies the calculations in the gluon 
sector. As always, physical parameters as Debye and 
Meissner masses do not depend on A. 
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III. INHOMOGENEOUS GLUON where Sfi = S^ c , even though m 2 M > in this interval. 

CONDENSATE AND INDUCED IN-MEDIUM Once Sfi > <5^ c , the Meissner mass square m 2 M becomes 

MAGNETIC FIELD negative itself, so the two terms of the effective Meissner 

mass give tachyonic contributions. Although the insta- 

To investigate the condensation phenomenon triggered bility in the region where m 2 M > is a direct consequence 

by the tachyonic modes of the charged gluons, we can re- of the pairing mismatch, it formally resembles the insta- 

strict our analysis to the gauge sector of the mean-field ef- bility produced by an external magnetic field on an other- 

fective action (P) that depends on the charged gluon fields wise stable color superconductor [lfj ■ The formal role of 

and rotated electromagnetic field. For a static solution, the magnetic field is played here by /x 8 . Borrowing from 

one only needs the leading contribution of the polariza- the experience gained in the case with external magnetic 

tion operators in the infrared limit (p Q = 0,\~]?\ — > 0). field, we would expect that this Meissner instability will 

Under these conditions, the gauge sector of the effective also be removed through the spontaneous generation of 

action can be written as an inhomogeneous gluon condensate (Ki) capable to in- 

/ ^ duce a rotated magnetic field, thanks to the anomalous 

d 4 x {--(J ll y)' 2 |n M if„ - IVyK^ 2 magnetic moment of the spin-1 charged particles. Having 

4 2 kkjg m mmclj we propose the following ansatz 

- [mlt S^iSui + (Ms - At3) 2 <V + nf^AK^Kl 



where m 2 M = 2a ^ [1 — ^-], is the Meissner mass with 

-p = fj, - + if and a s = £ Q, the last term in © 
comes from the T g term in fl}, taken in the 't Hooft 
gauge with gauge-fixing parameter A, II M = <9 M — iqAp, 

and ffj, u = d^Ay — <9„A M . In |(5J) the Debye mass tud 
was not included since it will be no substantial for our 
derivations. The chemical potential /13, although is zero 
in the gapped phase, should be in principle taken into 
account in the analysis of the new phase, since the K- 
condensate will break the remaining SU (2)c symmetry. 

At this point let us consider for a moment that we have 
an external rotated magnetic field H. In this case the ef- 
fective action @ becomes that of a spin-1 charged field 
in a magnetic field (see for instance [Hj|). If we work 
in the region where m 2 M — > 0, the ground state is 
stable = in the stable region) . Under these circum- 
stances, we know (Tpj that when qH > qH c = m 2 M — 
the effective magnetic mass of one of the charged gluon 
modes becomes imaginary due to the anomalous mag- 
netic moment term iqf^K^Kl- This field-induced in- 
stability triggers the formation of a gluon-vortex state 
characterized by the antiscrcening of the magnetic field. 

Now, let us go back to the situation of interest in the 
present work, that is, a system with no external magnetic 
field. As usual in theories with zero-component gauge- 
field condensates [l7| . us gives rise to a tachyonic mass 
contribution. This tachyonic contribution creates a re- 
gion of instability in a density interval for which m 2 M is 
still positive. Let us elaborate on this point. Using the 
results of [3] , one can easily realize that when the density 
is decreasing, coming from values where 8/j < A/ y2, the 
different chemical potentials accordingly change in such 
a way that the pairing mismatch 5fi increases. Conse- 
quently, the Meissner mass m 2 M decreases and, at some 
density value, it crosses the line of fig, leading to a tachy- 
onic effective Meissner mass {m 2 M — //§ < 0) for the range 
of densities between the crossing point and the value 



where we took advantage of the SU(2) C symmetry to 
write the (A'i)-doublet with only one nonzero component. 
Since in this ansatz the inhomogeneity of the gluon con- 
densate is taken in the (x, y)-plane, it follows that the cor- 
responding induced magnetic field will be aligned in the 
perpendicular direction, i.e. along the z-axes, (/12) = B. 
The part of the free energy density that depends on the 
gauge-field condensates, T g = T — T n o, with T n o = 
— r = J7 denoting the system free-energy density in the 
absence of the gauge-field condensate (G — 0,B = 0), 
is found, after fixing the gauge parameter to A = 1 and 
using the ansatz j3| in ((2]), to be 

T g = ^-2G*fl 2 G + 2g 2 \G\ 4 

-2[2qB + ^ 3 + ^) 2 + m 2 M }\G\ 2 (4) 

From the neutrality condition for the 3 r<i -color charge it 
is found that /13 = /is. The fact that /13 gets a finite value 
just after the critical point m 2 M — /ig = is an indication 
of a first-order phase transition, but since fis is para- 
metrically suppressed in the gapped phase by the quark 
chemical potential /is ~ A 2 //z [J], it will be a weakly 
first-order phase transition. Henceforth, we will consider 
that [13 = fj,s in an d work close to the transition 
point 5/i > 5/j, c which is the point where m 2 M continu- 
ously changes sign to a negative value. For very small, 
negative values of ra 2 M , the gluon condensate and the in- 
duced magnetic field should be very small too, thereby 
facilitating the calculations. 

Minimizing (0J with respect to G gives 

-U 2 G-(2qB + m 2 M )G + 2g 2 \G\ 2 G = (5) 

Eq. (O is a highly non-linear differential equation that 
can be exactly solved only by numerical methods. Never- 
theless, we can take advantage of working near the transi- 
tion point, where we can manage to find an approximated 
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solution that will lead to a qualitative understanding of 
the new condensate phase. With this aim, and guided by 
the experience with external fields, where the solution 
is always such that the kinetic term — Ii u K il \ 2 

is approximately zero near the transition point, we will 
consider that when Sfj, ~ <5/i c our solution will satisfy 
the same condition. Hence, we will look for a minimum 
solution satisfying 

n 2 G + qBG~0. (6) 

With the help of ([6]) one can show that the minimum 
equation for the induced field B takes the form 

2q\G\ 2 - B ~ (7) 

The relative sign between the two terms in Eq. ((TJ) im- 
plies that for \G\ ^ a magnetic field B is induced. 
The origin of that possibility can be traced back to the 
anomalous magnetic moment term in the action of the 
charged gluons. This effect has the same physical root 
as the paramagnetism found in Ref. [101 ]; where con- 
trary to what occurs in conventional superconductivity, 
the resultant in-medium field B becomes stronger than 
the applied field H that triggers the instability. The an- 
tiscreening of a magnetic field by the condensation of 
charged spin-1 fields has been also found in the context 
of the electroweak theory [13] ■ 

Using the minimum equations ([5]) and ([7]) in ((4]), we 
obtain the condensation free-energy density 

T g ^-2(g 2 -t)\G\" (8) 

The hierarchy between the strong (a) and the electro- 
magnetic (q) couplings implies that T c < 0. Therefore, 
although the induction of a magnetic field always costs 
energy (as can be seen from the positive first term in f4|), 
the field interaction with the gluon anomalous magnetic 
moment, produces a sufficiently large negative contribu- 
tion to compensate for that increase. Consequently, as 
seen from (jHJ , the net effect of the proposed condensates 
is to decrease the system free-energy density. 

It follows from Eqs.(|5])-((7j) that near the phase transi- 
tion point the inhomogeneity of the condensate solution 
should be a small but nonzero correction to a leading 
constant term 

|G| 2 ~ A^il^q 2 + G(m 4 M )/(x, y), (9) 

qB ~ A g/ gm 2 M + 0{m 4 M )g{x, y). (10) 

with A fl/? = (sV^-l)- 1 . 

The explicit form of the inhomogeneity can be found 
from ©, which can be written in polar coordinates as 




In the above equation we approached the rotated mag- 
netic field by its leading in (TTT)]) , used the symmetric 
gauge Ai = —(B/2)eijXj, and introduced the characteris- 
tic length £ 2 = 1/ ' k g /qm 2 M . Using just the leading contri- 
bution of B in pip is a consistent approximation if we si- 

2 

multaneously drop the term and restrict the solution 
to the domain r<^. Notice that this domain is indeed 
a large region because near the critical point the magni- 
tude of rriM is very small, hence £ becomes very large. 
The most symmetric solution of Eq. (jTTJ) is the one that 
preserves the SO (2) symmetry in the (x, y) plane. Hence, 
proposing a solution of the form G(r, 9) ~ R(r)e tx , with 
X a constant phase, the equation for R(r) reduces to 
\rd r (rd r ) + |j-]i?(r) = 0. It is solved by the Bessel func- 
tion of the first kind Jo (»"/£)■ Then, the gluon condensate 
can be written as G(r) — (l/y/2q£)Jo(r/£) expi\, which 
is consistent with (|9|) , as in the domain of validity of this 
solution (r <C £) the series can be approximated by its 
first terms. Accordingly, the modulus of the condensate 
square satisfies 

— 2 1 r2 

The improved solution for B is found substituting (TT2"|) 
back into . The induced field B is homogeneous in the 
z-direction and inhomogeneous in the (x, y)-plane. 

As always, the choosing of a particular gauge condi- 
tion (A = 1 in this case) is dictated by convenience. As 
occurs in other gauge theories, the particular form of the 
condensate solution G(r, 9) and the effective action (or 
the effective potential in the case of constant conden- 
sates) are gauge-dependent. There is nothing wrong in 
that, as none of these quantities is physically measur- 
able. What is gauge independent is the effective action 
evaluated in the minimum solution. Therefore, if the 
result (fT2|) is plugged back into (jSJ) , one obtains a gauge- 
independent quantity representing the free-energy den- 
sity of the ground state. 

One may wonder whether this inhomogeneous gluon 
condensate forms a vortex state. To answer this question 
we can compare our results with the case with external 
magnetic field [Toj ] . For this we should notice that the 
mathematical problem we have just solved is formally 
similar to that where the instability is induced by a weak 
external field. This would be the situation when the 2SC 
system approaches the transition point from the stable 
side (real effective magnetic mass) and the external mag- 
netic field is slightly larger than a very small but still 
positive effective mass square H > H c = m 2 M — /x| <C 1. 
We know that at large m 2 M —n\ the condensate solution is 
a crystalline array of vortex-cclls with cell's size ~ £ <C 1 . 
At smaller m 2 M — /xf the lattice structure should remain, 
but with a larger separation between cells, since in this 
case £ 3> 1. However, the use of a linear approximation 
to solve the equations in this case only allows to explore 
the solution inside an individual cell (r <^ £). This is 
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the same limitation we have in the linear approach fol- 
lowed in the present paper. Therefore, we expect that 
when the nonlinear equations will be solved, the vortex 
arrangement will be explicitly manifested. 

The idea of removing the 2SC chromomagnetic in- 
stability with the help of a gluon condensate has been 
previously considered in Ref. Q. In that work a ho- 
mogeneous gluon condensate with no rotated magnetic 
field was proposed. Near the transition point the homo- 
geneous condensate produces a contribution to the con- 
densation energy density that can be approximated as 
~T g ~ ~^j(m 2 M ) 3 /m 2 , with m 2 = 4a 2 s JI 2 /3ir. In this 
case, the condensation energy in a circle of radius c£, 
with c being a dimcnsionless constant satisfying c < 1, is 

~F ~ 2%. On the other hand, for the G — B 

condensate considered in the present paper, the corre- 
sponding condensation energy is F g ~ — c 2 nm 2 iI /2q 2 , 
which is found integrating the free-energy density ([8]) 
evaluated in the solution (fT2|) in the circle of radius 
c£. Taking into account that near the transition point 
m 2 M <C m 2 [4| , one can see that the inhomogeneous gluon 
condensate is energetically favored over the homogeneous 
one. 

We should stress that as 5fi — Sfj, c becomes much larger, 
the contribution of G and B in the quasiparticle propaga- 
tor S cannot be neglected, thereby affecting the gap equa- 
tion (as shown in [18j , a sufficiently strong magnetic field 
can qualitatively influence the gap). As a consequence, 
the inhomogeneity of the gluon condensate will be nat- 
urally transferred to the diquark condensate. Hence, a 
LOFF-type phase may appear as a back reaction of the 
inhomogeneous gluon condensate on the gap solution. It 
is natural to expect that the additional reduction in free- 
energy due to the G — B condensation will make this 
phase energetically favored over a pure LOFF one 0] ■ 

We anticipate that a G — B condensate will likely 
remove the chromomagnetic instability in gCFL too. 
In the gCFL case there are four gluons with tachyonic 
masses. Following the results of the last paper in Ref. 

the four tachyonic gluons are A\, A% and two com- 
binations of A3, As and A 1 , A third combination of A3, 
As and A 1 is massless, hence it represents the rotated 
electromagnetic field in the gCFL phase. Notice that, in 
the unstable gCFL phase the gluons A\ and A2 acquire 
rotated charge since they couple to the new rotated elec- 
tromagnetic field through its A3 component. This im- 
plies that Ai and A2 are analogous to the charged glu- 
ons that become tachyonic in the 2SC case. If A\ and 
A2 condense in an inhomogeneous condensate following 
the same mechanism we found in the present paper, this 
condensate could induce a rotated magnetic field and also 
give real masses to the two combinations of A3, As and 
A 1 that were tachyonic. Of course, it is possible too 



that in addition to the condensation of the charged glu- 
ons, the other two tachyonic combinations also condense 
in the gCFL case. The only way to find out which sce- 
nario is the most energetically favored is by doing explicit 
calculations. However, exploring this idea will be more 
challenging, as one will have to deal with a strongly first- 
order phase transition. 



IV. ASTROPHYSICAL CONNOTATIONS 

A common characteristic of neutron stars is their 
strong magnetization. Their surface magnetic fields 
range from B = 1.7 x 10 8 G (PSR B1957+20) up to 
2.1 x 10 13 G (PSR B0154+61), with a typical value of 
10 12 G [19(. There are observational discoveries of even 
strongly magnetized stellar objects- known as magnetars- 
with surface magnetic fields of order B ~ 10 14 — 10 16 G 
[20l |. In the core of these compact objects, the field may 
be considerably larger due to flux conservation during 
the core collapse. By applying the scalar virial theorem 
it can be shown that the interior field can reach values of 
order B ~ 10 18 G [p. 

The observed stellar magnetic fields are supposed to be 
created by a magnetohydrodynamic dynamo mechanism 
that amplifies a seed magnetic field due to a rapidly rotat- 
ing protoneutron star. Thus, the standard explanation of 
the origin of the magnetar's large magnetic field implies 
that the rotation should have a spin period < 3ms. Nev- 
ertheless, this mechanism cannot explain all the features 
of the supernova remnants surrounding these objects [22I ]. 

It has been found [4| that when Sfj, — A in the neu- 
tral 2SC model, the absolute value of the magnetic mass 
becomes of order m g . This implies that the inhomoge- 
neous gluon condensate found in our paper could pro- 
duce a magnetic field of order of magnitude in the range 
- 10 16 - 10 17 G. As discussed in Refs. the pos- 

sibility of generating a magnetic field of such a large 
magnitude in the core of a compact star, without re- 
lying on a magneto- hydrodynamics effect, can be an in- 
teresting alternative to address the main criticism [22| to 
the observational conundrum of the standard magnetar's 
paradigm [2fl |. On the other hand, to have a mechanism 
that associates the existence of high magnetic fields to CS 
at moderate densities can serve to single out the magne- 
tars as the most probable astronomical objects for the 
realization of this high-dense state of matter. 
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